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ABSTRACT 

We  report  on  the  optical  properties  of  two  dimensional  (2D)  photonic  crystals  (PCs)  deeply 
etched  in  an  InP/GalnAsP  step-index  waveguide.  Transmission  (T)  measurements  through 
simple  PC  slabs  and  through  one-dimensional  (ID)  Fabry-Perot  (FP)  cavities  between  PC 
mirrors  are  reported  and  compared  to  theory.  A 2D  finite  difference  time-domain  (FDTD) 
method  combined  to  a phenomenological  out-of-plane  loss  model  is  used  to  assess  different  loss 
contributions.  The  PC  optical  properties  are  deduced  from  the  FP  peak  analysis.  The  origin  of  the 
high  T level  observed  inside  the  stopgap  is  investigated. 


INTRODUCTION 

In  photonic  crystal  (PC)  structures  the  periodic  arrangement  of  different  elements  with 
strong  dielectric  contrast  affects  the  properties  of  photons  forbidding  light  propagation  along 
specific  directions  within  relatively  large  energy  bands  known  as  photonic  hand-gaps  (PBGs) 
[1,2].  In  the  last  few  years  it  has  been  proved  that,  although  an  omnidirectional  PBG  is  only 
possible  in  a three-dimensional  (3D)  PC,  a two-dimensional  (2D)  PC  combined  to  a vertical  step- 
index  waveguide  offers  enough  light  control  for  integrated  optics  applications  [3].  The  potentials 
of  this  approach  have  been  successfully  demonstrated  both  in  GaAs-based  structures  (k~  1.0 
flm)  [3]  and  in  InP-based  systems  at  telecommunication  application  wavelengths  (k  ~ 1.5  pm) 
[4].  In  this  paper  we  describe  the  optical  characterization  of  quasi-2V  PCs  deeply  etched  into  an 
InP/GalnAsP  slab  waveguide.  The  internal  light  source  (ILS)  technique  is  used  to  measure  the 
transmission  (T)  spectra  through  simple  PC  slabs  and  one-dimensional  (ID)  Fabry-Perot  (FP) 
cavities  formed  between  two  PC  mirrors.  A 2D  finite  difference  time-domain  (FDTD)  method  is 
used  to  best  fit  the  experimental  data.  The  exact  amount  of  out-of-plane  losses  is  assessed  and 
the  hole  depth/shape  contributions  are  analyzed.  The  optical  properties  of  the  PC  structures  are 
determined  and  the  origin  of  the  high  T level  detected  inside  the  stopgap  is  discussed. 


EXPERIMENT 

Quasi- 2D  PCs  consisting  of  a triangular  lattice  of  air  holes  were  fabricated  on  nominally 
undoped  GalnAsP/InP  heterostructures  grown  by  metal  organic  vapor  phase  epitaxy  (MOVPE) 
on  n-InP  substrates  (see  figure  1(a)).  Moderate  filling  factor  values  (f~  0.30)  were  chosen  in 
order  to  minimize  out-of-plane  losses  [3,4],  Such  low/values  result  in  a full  gap  only  in  the 
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transverse  electric  (TE)  polarization  direction.  The  PC  structures  were  etched  normally  to  an 
InP/GalnAsP/InP  slab  waveguide  containing  a built-in  light  source.  The  waveguide  core  is  a 434 
nm-thick  Gan.24Ino.76Aso.52Pn.4x  layer  lattice  matched  to  InP  and  sandwiched  between  a 200  nm- 
thick  InP  cap  layer  and  a 600  nm-thick  InP  buffer  layer.  In  the  spectral  region  of  interest  (/.<?.,  X 
~ 1 550  nm)  the  values  = 3.35  and  «,/,„/  = 3. 17  arc  assumed  for  the  refraction  index  of 
GalnAsP  and  InP,  respectively.  The  resulting  structure  is  a multimode  waveguide  with  a 
fundamental  guided  TE  mode  corresponding  to  an  effective  index  «,#  = 3.24  [4],  The  light  source 
consists  of  two  different  GalnAsP  strain-compensated  quantum  well  (QW)  packages  separated 
by  a 30  nm  spacer  and  located  approximately  in  the  middle  of  the  guiding  structure  between  two 
181  nm-thick  barrier  layers.  Due  to  the  superimposition  of  the  two  QW  signals  (X.|  = 1.47  pm 
and  X2  = 1.55  pm)  a 100  nm-wide  PL  peak  centered  at  1500  nm  is  obtained  [4]. 

Electron  beam  (e-beam)  lithography  was  used  to  define  the  PC  patterns,  while  dry  etching 
techniques  were  applied  to  transfer  the  pattern  into  the  semiconductor  heterostructure.  In  order  to 
form  a durable  SiOi  etch  mask  for  the  subsequent  etch  process  the  procedure  described  in  [4,5] 
was  adopted.  Ar/CL  chemically  assisted  ion  beam  etching  (CAIBE)  [4]  and  Ch/Ar-based 
electron  cyclotron  resonance-RIE  (ECR-RIE)  [5]  were  used  to  pattern  sample  #1  and  sample  #2, 
respectively.  As  sketched  in  figure  1(a),  the  obtained  hole  shape  is  conical  with  nearly  vertical 
walls  close  to  the  surface  and  a strongly  tapered  bottom.  Figures  1(b)  and  (c)  show  the  cross- 
section  view  scanning  electron  microscopy  (SEM)  micrographs  of  two  PC  test  structures  for 
sample  #1  and  sample  #2,  respectively.  In  sample  #1  the  holes  have  oblique  walls  already  inside 
the  waveguide  core  (-2)  and  a slight  bending  towards  the  hole  bottom  is  observed.  Nevertheless, 
the  hole  pattern  is  "ordered":  no  dispersion  in  the  hole  depth  : = zi  + A:  is  found  when  the  hole 
diameter  (D)  is  kept  constant,  r values  of  2.55  pm  are  obtained  for  the  larger  holes  [4],  On  the 
contrary,  in  sample  #2  fluctuations  appear  in  the  hole  depth  yielding  <->  > 3 pm,  while  no 
bending  is  present  and  the  holes  have  straight  walls  inside  -2. 

The  ILS  technique  described  in  [4]  was  used  to  characterize  the  PC  optical  properties. 


Figure  1.  (a)  Sketch  of  the  InP/GalnAsP  waveguide  (WG)  hctcrostructure.  The  hole  shape 
model  used  to  calculate  the  loss  parameter  value  e"  is  illustrated.  z2  is  the  thickness  of  the  cap  + 
WG  layer  where  the  holes  are  assumed  to  be  cylindrical,  while  Az  is  the  cone  depth.  The  cone 
half  angle  a is  also  evidenced,  (b)  and  (c):  edge  view  SEM  micrographs  of  test  photonic  crystal 
(PC)  structures  for  sample  #1  and  sample  #2,  respectively.  The  images  were  taken  before  the 
removal  of  the  SiCF  mask.  In  the  inset  of  figure  (c)  the  cylindrical  shape  of  sample  #2  holes 
inside  Z2  is  evidenced. 
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The  PL  excited  inside  the  GalnAsP  QWs  is  used  as  a built-in  probe  beam.  Part  of  the  PL  signal 
propagates  parallel  to  the  surface  as  a guided  mode  and  interacts  with  the  PC  structure.  The 
absolute  PC  T spectrum  is  obtained  normalizing  the  spectrum  detected  after  transmission  through 
the  PC  structure  with  respect  to  the  spectrum  collected  in  a non-patterned  region  of  the  sample. 
Details  on  the  experimental  set-up  and  on  the  measurement  schemes  adopted  for  simple  PC  slabs 
and  ID  FP  cavities  are  given  in  [4],  Here  we  remind  that,  since  ILS  measurements  on  a single  PC 
structure  yield  the  T spectrum  only  in  a narrow  interval  (~  100  nm),  the  "lithographic  tuning" 
approach  has  to  be  applied  in  order  to  determine  the  PC  optical  properties  within  a wide  spectral 
range.  As  the  QW  emission  wavelength  is  fixed,  the  scaling  property  of  PCs  f 1]  is  used  and  PC 
structures  with  different  periods  a are  measured.  If/  is  kept  constant  the  whole  PBG  can  be 
explored  as  a function  of  the  reduced  frequency  u = a /X  [3,4]. 


RESULTS  AND  DISCUSSION 


T spectra  through  PC  slabs  along  both  TM  and  UK  orientations  and  for  TE  polarization  are 
shown  in  figures  2 (a)  and  (b)  for  sample  #1  (10  rows  thick  slab)  and  #2  (8  rows),  respectively. 
The  layout  of  simple  PC  slabs  is  shown  in  figure  2 (c).  Both  TM  and  TK  oriented  structures  were 
fabricated  with  periods  a ranging  from  240  to  500  nm  for  sample  #1  (A a = 20  nm)  and  from  240 
to  640  nm  for  sample  #2  (A a = 20  nm).  Well-defined  PBGs  appear  in  each  spectrum  with  steep 
dielectric  band  edges  at  u = 0. 19  and  0.22  for  TM  and  TK  orientations,  respectively.  On  the  other 
hand,  due  to  the  strong  influence  of  out-of-plane  losses  on  T at  the  air  band  edge  [3],  smooth  T 
bands  appear  at  u ~ 0.28  and  0.30  for  FM  and  UK  orientations,  respectively.  T values  between 
80%  and  90%  are  reached  for  u < 0.20,  while  T < 40%  for  u > 0.30.  Finally,  an  average  T = 3 - 
5%  is  observed  inside  the  PBG  for  both  samples. 

The  experimental  T spectra  are  compared  to  theoretical  spectra  calculated  with  a 2D  FDTD 
model  [6].  The  effect  of  the  vertical  confinement  on  the  light  propagating  through  the  PC  slabs 
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Figure  2.  TE  transmission  spectra  through  TM  and  TK  oriented  photonic  crystal  (PC)  slabs 
for  (a)  sample  #1  (10  rows)  and  (b)  sample  #2  (8  rows),  respectively.  Experimental  spectra 
(black  lines)  are  compared  with  2D  finite  difference  time-domain  (FDTD)  calculated  spectra 
(gray  lines),  (c)  Sketch  of  the  typical  layout  of  simple  PC  structures  along  TM  and 
TK  orientations.  Each  slab  is  characterized  by  the  period  a (or  hole  diameter  D)  value.  Hairlines 
show  the  single  "atomic"  plane  for  both  orientations. 
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was  taken  into  account  assuming  n = n,g  = 3.24  for  the  dielectric  matrix.  Out-of-plane  losses 
were  cast  into  the  2D  calculation  by  introducing  the  imaginary  e”  as  a phenomenological  loss 
parameter  [4,7]./ and  e"  were  chosen  as  free  parameters  of  the  fit.  The  best-fitted  curves  are 
reported  in  figures  2(a)  and  (b),  while  in  table  I the  obtained /and  e"  values  are  listed.  While  for 
sample  #2  both  TM  and  TK  spectra  arc  fitted  assuming  the  same  /value,  for  sample  #1  / (TK)  > / 
(TM).  The  latter  discrepancy  is  due  to  the  hole  squeezing  along  the  TK  direction  originating  in 
the  e-beam  lithography  of  sample  #1  [4].  Different  e"  values  have  to  be  used  to  fit  cither  the 
dielectric  or  the  air  T band  edge  in  both  samples.  As  shown  in  [4],  e"  <=  (iC  xf):  if/ is  constant 
while  increasing  «,  higher  e"  values  have  to  be  introduced.  Moreover,  when  it  is  fixed,/ (rK)  >/ 
(rM)  implies  £"  (rK)  > e"  (TM).  Finally,  the  sample  #2  e"  values  arc  always  higher  than  the 
corresponding  values  of  sample  #1 . The  different  out-of-plane  loss  level  in  the  two  samples  has 
to  be  ascribed  to  differences  in  the  hole  morphology  and  in  the  PC  pattern  disorder  (see  below). 

The  e"  values  obtained  from  the  best  fit  of  the  air  band  edge  were  analyzed  in  the  framework 
of  a loss  model  where  the  total  loss  parameter  is  expressed  as  £"  = £",,„  + £"/„,/,.  [8,9],  £"/,„/<. 
contains  contributions  from  the  finite  etch  depth  and  from  the  conical  hole  shape,  while  £",-„, 
accounts  for  the  intrinsic  losses  corresponding  to  the  ideal  case  of  infinitely  deep  holes.  Since  for 
these  InP-based  PC  structures  e",„,  = 0.01  - 0.02  [4],  the  values  (#1 ) = 0. 105  ± 0.045  and 
(#2)  =0.305  + 0.005  derive  for  sample  #1  and  #2,  respectively.  Following  the  procedure 
illustrated  in  [9]  and  adopting  the  model  hole  shape  sketched  in  figure  1(a).  the  theoretical  value 
can  be  calculated.  If  the  exponential  decaying  field  profile  !/:)  = A exp  (Kz)  is  assumed  in 
the  substrate,  £";,„/<•  is  a function  of  K and  of  the  cone  slope  a in  the  first  L,/<r„v  = AT1  [9],  From  the 
SEM  micrograph  analysis  a ~ 2.5°  ± 0.5°  for  both  samples  and  £"/,„/,,  = 0.30  is  obtained.  While 
the  latter  value  perfectly  agrees  with  (#2),  it  is  2 - 4 times  higher  than  (#1 ).  As 
mentioned  in  [4],  the  discrepancy  between  the  calculated  and  the  experimental  £"/,„/,.  (#1)  values 
originates  from  the  complex  hole  geometrical  shape  (e.g.,  hole  bending,  oblique  walls  inside  the 
core)  which  does  not  fit  the  assumed  model  shape.  On  the  contrary,  when  holes  keep  straight 
walls  inside  the  cap  + waveguide  layer  (rj),  the  different  loss  contributions  can  be  deduced  from 
the  e"  experimental  values  and  £"/,„/,■  can  be  easily  modeled.  Extrapolating  from  sample  #2  data, 
two  perspectives  emerge  which  can  be  considered  realistic  for  the  actual  etching  techniques. 

First,  it  is  found  that  the  minimum  loss  limit  £"„,,  could  be  reached  for  a < 0.5°,  i.e.  for  holes  with 
almost  straight  walls  in  some  s inside  the  substrate.  On  the  other  hand,  a less  critical 
reduction  in  the  cone  slope  would  make  possible  to  get  a sufficiently  low  loss  level  for  good 
performance  PC  applications  (e.g.,  t"  <0.1  for  a < 1 .5°)  [8], 

In  order  to  investigate  the  PC  optical  properties  inside  the  PBG,  ID  FP  cavities  were 
fabricated  for  both  samples  using  4 rows  thick  TM-orientcd  PC-based  mirrors  (a  = 300  - 480  nm 


Table  I.  Photonic  crystal  (PC)  slabs  : / and  £"  values  (FDTD  fit).  ID  Fabry-Perot  (FP)  cavities  : 
optical  properties  of  a single  4 row  PC  mirror  and  cavity  parameter  values  (Airy's  function  fit). 


Sample 

PC  slabs  - FDTD  fit 

ID  - FP  Cavities  - Airy's  function  fit 

N 

/ 

E"  (Diel.  B.) 

E"  (Air  B.) 

T (%) 

R(%) 



L (%)  Lp la 

m 

Q 

#1 

10 

TM 

0.30 

0.04 

0.08 

a = 340  nm  ; Wla  = 1 .9 

TK 

0.32 

0.06 

0.16 

14 

64 

22  1 0.36 

4 

28 

TM 

0.30 

0.08 

0.32 

a = 340  nm  ; Wla  = 2.0 

#2 

8 

FK 

0.30 

0.08 

0.32 

10 

80 

10  | 0.40 

4 

56 
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and  Aa  = 20  nm)  separated  by  a spacer  with  a normalized  width  W I a = 1.7  - 2.2  [4].  A typical  T 
spectrum  through  a ID  FP  cavity  is  shown  in  figure  3(a)  for  sample  #2  (W / a = 2.0).  The  cavity 
layout  is  sketched  in  figure  3(b).  The  two  peaks  appearing  in  the  T spectrum  at  u = 0.22  and  u = 
0.29  are  due  to  first  and  second-order  FP  resonances,  respectively.  The  experimental  spectrum  is 
compared  to  the  calculated  2D-FDTD  T spectrum.  The  TM  values  derived  from  the  best  fit  of 
simple  PC  slab  air  band  edges  were  adopted  for/ and  e",  while  the  value  W / a = 1 .97  was 
assumed  for  the  cavity  width.  As  for  simple  PC  slabs,  FDTD  calculated  spectra  do  not  agree  with 
the  experimental  data  inside  the  PBG:  the  theoretical  T value  inside  the  gap  (T  = 10  2 - 104  %)  is 
higher  than  the  detected  value  (T  ~ 3 - 5 %).  While  the  experimental  FP  peaks  are  located  at 
energies  consistent  with  the  calculated  values,  they  show  higher  peak  T and  lower  quality  factors 
Q than  the  simulated  peaks.  As  shown  in  the  inset  of  figure  3(a)  the  Airy's  formula  for  planar 
resonators  is  used  to  fit  the  first-order  FP  resonances  [4]  yielding  the  optical  properties  of  a 
single  4 row  PC  mirror  (see  table  I). 

The  comparison  between  the  separate  analysis  of  simple  PC  slab  and  1 D FP  cavity  spectra 
for  both  samples  completes  the  picture  given  in  [4],  While  the  fit  of  the  T air  band  edges  yields 
e"(#2)  > e"(#l),  from  the  best  fit  of  FP  resonances  inside  the  PBG  the  total  loss  value  L (#2)  is 
found  to  be  lower  than  the  corresponding  L (#1)  value.  Moreover,  in  sample  #2  the  T level  inside 
the  stopgap  decreases  while  reflectivity  (R)  increases  to  80%  giving  Q (#2)  = 2 xQ  (#1).  In 
order  to  explain  this  apparent  discrepancy  we  remind  that  two  different  effects  may  affect  light 
propagation  through  a PC  structure  [4],  First,  the  finite  hole  depth  contributes  to  increase  out-of- 
plane  scattering.  On  the  other  hand,  when  entering  the  PC,  the  hole  pattern  is  felt  as  a low 
refractive  index  contrast  region  and  light  is  pushed  below  the  guiding  layer  ( modal  conversion). 
When  energies  outside  the  stopgap  are  considered,  light  propagation  is  allowed  and  Bloch  modes 
propagate  into  the  PC.  As  pointed  out  in  [10],  Bloch  modes  are  strongly  influenced  by  defects 
and/or  disorder  where  out-of-plane  scattering  can  occur.  Therefore,  due  to  dispersion  in  the  cone 
depth  Az  (see  figure  1(c))  which  contributes  to  disorder,  out-of-plane  losses  are  found  to  be 
higher  in  sample  #2  than  in  sample  #1.  On  the  contrary,  the  second  effect  becomes  predominant 
inside  the  stopgap.  Light  propagation  is  forbidden  and  "evanescent"  modes  penetrate  the  PC. 
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Figure  3.  (a)  Transmission  spectrum  through  a ID  Fabry-Perot  (FP)  cavity  between  two  TM 
oriented  4 rows  photonic  crystal  (PC)  mirrors  separated  by  a spacer  W.  Arrows  indicate 
transmission  peaks  related  to  first  and  second  order  FP  resonances.  Experimental  spectra  (black 
lines)  are  compared  with  2D  finite  difference  time-domain  (FDTD)  calculated  spectra  (gray 
lines).  The  Airy's  function  best-fit  (gray  dashed  line)  of  the  first  order  peak  for  a - 340  nm  is  the 
shown  in  the  inset,  (b)  Sketch  of  the  typical  layout  for  ID  FP  cavity.  The  actual  cavity  width  W ' 
taking  into  account  the  finite  penetration  depth  Lp  of  the  field  inside  the  mirrors  is  shown. 
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While  these  modes  are  insensitive  to  defects  and  disorder,  they  arc  strongly  influenced  by  the 
modal  conversion  and  light  propagation  is  affected  by  the  single  hole  shape  parameters  (e.g.,  the 
hole  depth).  Thus,  since  <c>  (#2)  > z (#!),  inside  the  stopgap  light  transmission  is  lowered  in 
sample  #2  with  respect  to  sample  #1 , while  R and  cavity  Q values  result  to  be  improved. 


CONCLUSIONS 

We  have  reported  on  the  optical  properties  of  2D  PCs  etched  in  InP/GalnAsP  slab 
waveguide  structures.  The  combination  between  structural  analysis  and  optical  characterization 
by  ILS  experiments  makes  possible  to  evaluate  which  improvements  are  necessary  in  the 
fabrication  process.  Two  different  regimes  are  identified  for  light  propagation  inside  the  PC 
structure.  When  light  with  energies  outside  the  PBG  is  considered,  out-of-plane  losses  are  the 
limiting  factor  for  PC  performances.  Improvements  in  the  hole  shape  and  in  the  ordering  of  the 
PC  pattern  should  allow  one  to  lower  the  scattering  towards  the  substrate.  On  the  other  hand, 
when  entering  the  stopgap,  the  hole  "physical"  parameters  (e.g.,  the  hole  depth)  strongly  affect 
the  PC  optical  properties. 
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